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ABSTRACT 


Force  measurements  were  made  to  investigate  the  effect  of  wall 
temperature  on  sphere  drag  in  low-density  supersonic  flow.  Drag  data 
were  obtained  by  fixing  flow'  conditions  and  changing  only  the  sphere 
wall  temperature.  Conditions  covered  in  the  experiment  include 
3. '3  k  Mjb  <  3.  7,  4.  0  <  Re^  <  10.  0,  0.  5  <  Kn^  <  1.  3,  and  1.2  <  Tw/T00  <  3.0. 
A 'means  of  correlating  drag  data  in  transition  flow  is  suggested,  and  its 
validity  assessed  by  applying  the  correlation  parameter  to  this  and 
several  other  sets  of  experimental  data. 
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NOMENCLATURE 

A 

Cross-sectional  area  of  sphere,  nr2 

CD 

Drag  coefficient,  D /( 1  /  2  pa>Uca2A) 

CP 

Specific  heat  of  gas  at  constant  pressure 

Cv 

Specific  heat  of  gas  at  constant  volume 

D 

Drag 

d 

Sphere  diameter 

erf(x) 

x  2 

Error  function,  (2/J~n)  f  e_t  dt 

o 

f 

Velocity  distribution  function 

g 

Relative  speed  of  molecules 

H 

Local  total  enthalpy 

h 

Defined  by  Eq.  (5e) 

I 

Defined  by  Eq.  (8) 

Kn 

Knudsen  number  based  on  sphere  diameter,  X/d 

Kw,  d 

Kr^KyJV  Sw  +  1) 
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M 

m 

N 

n 

P 

PP 

PP 

q 

R 

Re 

Reo,  r* 
Re* 

r 

r* 

S 

sr 

Sw 

S. 

T 

T 

x  r 

t 
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u 

V 

V 

w 

J3 

7 


Mach  number 
Mass  of  a  molecule 
Collision  frequency 
Number  density  of  molecules 
Pressure 

Total  pressure  immediately  downstream  of  a  normal 
shock 

Measured  pitot  pressure 

Dynamic  pressure,  pU^/2 
Gas  constant 

Reynolds  number  based  on  sphere  diameter 
Nozzle  reservoir  Reynolds  number,  pQ  2HQ  r */p0 
p  U  d/jj 

1  CD  CD  CD 

Radius  of  sphere 
Nozzle  throat  radius 
Speed  ratio,  UA2RT)1/2 
U00/(2RTr)1/2 

Ucd/(2RTw)^^2 

U<D/(2RTa))l/2 

Temperature 

Reference  temperature,  (TQ~  4  T+)/2 

Time 

Velocity 

Molecular  velocity  component 
Molecular  velocity  component 
Mean  molecular  velocity 
Molecular  velocity  component 
Angle  defined  in  Fig.  5 
Ratio  of  specific  heats,  C  /Cv 
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£  Defined  by  Eq.  (5c) 

i?  Defined  by  Eq.  (5b) 

X  Mean  free  path 

p  Viscosity 

5  Defined  by  Eq.  (5a) 

p  Gas  mass  density 

o  Effective  molecular  diameter,  (T2/T1)1/2  or 

SUBSCRIPTS 

1  Corresponding  to  class  1  molecules 

2  Corresponding  to  class  2  molecules 

d  Corresponding  to  sphere  diameter 

fm  Free-molecular  value 

o  Reservoir  (total)  condition 

w  Sphere  wall  condition 

*  Free -stream  condition 


SUPERSCRIPTS 


+ 


Denotes  incoming  molecules 
Denotes  outgoing  molecules 
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SECTION  I 
INTRODUCTION 


An  extensive  experimental  sphere  drag  investigation  has  recently 
been  completed  by  Bailey  and  Hiatt  (Ref.  1)  in  an  aeroballistic  range  of 
the  von  Karman  Gas  Dynamic  Facility  (VKF)  at  AEDC.  This  test  was 
conducted  to  obtain  sphere  drag  coefficients  for  data  reduction  of  the 
falling  balloon  experiments  (Ref.  2).  Their  measurements  covered  a 
Mach  number  range  of  0.  12  <  <  6.  39  and  Reynolds  number  range  of 

15  <  Re0  <  50,  300  at  a  T-w/Ta,  ratio  of  about  one.  Sphere  drag  measure¬ 
ments  were  also  made  in  the  Aerospace  Research  Chamber  (10V)  of  VKF 
as  a  part  of  this  program.  The  10V  data  are  reported  herein,  and  they 
cover  the  flow  conditions  3.  31  <  M,,,  <  3.  63,  4.  17  <  Re,,  <  9.  27, 

0.  55  <  Kn^  <  1.  29,  and  1.  20  <  Tw/TaD  <  2.  85.  In  addition,  the  nominal 
Mach  seven  (M7)  balance  data  reported  in  Ref.  3  were  reduced  using 
the  present  data  reduction  computer  program.  These  data  were  pre¬ 
viously  reduced  by  hand  calculations.  The  results  are  changed  less 
than  ±2  percent.  The  flow  conditions  of  the  M7  data  are  7.21  <  <  7.77, 

2.  85  <  Rem  <  8.  32,  1.  39  <  Kn^  <  3.  77,  and  3.  90  <  Tw/T„,  <  5.  44. 

The  primary  objective  of  the  10V  investigation  was  to  study  the 
effect  of  wall  temperature  on  sphere  drag  for  Knudsen  numbers  of 
order  one,  0(1).  A  difficulty  in  investigating  the  effect  of  Tw  is  that  the 
effect  of  other  flow  variables  must  be  either  known  or  eliminated.  Un¬ 
fortunately,  the  flow  conditions  required  to  produce  Kn^  -  0(  1)  for 
practical  model  sizes  (in  this  case  d  =  0.  25  in. )  are  difficult  to  calibrate. 
This  is  due  to  the  orifice  effect,  thermal  transpiration  effect,  and 
viscous  correction  on  pitot  pressure  measurements,  as  well  as  prob¬ 
lems  associated  with  calibrating  a  device  to  measure  pressures  near 
10  n  Hg. 

In  addition  to  flow  calibration  a  certain  amount  of  difficulty  is  in¬ 
volved  in  measuring  small  drag  forces  (in  this  case  the  drag  was  on  the 
order  of  100  mg).  The  tare  force  on  the  balance  sting  must  also  be 
measured  and  subtracted  from  the  total  force  measurement.  The  tare 
force  was  typically  from  30  to  35  percent  of  the  total  drag  force.  Force 
measurements  of  this  magnitude  are  sensitive  to  vibrations,  balance 
alignment,  and  thermal  environment,  any  one  of  which  can  cause 
meaningless  drag  measurements  if  not  controlled.  Obviously,  inter¬ 
ference  and  mechanical  problems  are  removed  when  free-flight  tech¬ 
niques  are  used;  however,  error  is  introduced  in  free-flight  data  by 
inaccuracies  in  the  sphere  position  measurements.  The  probable  magni¬ 
tude  of  error  and  scatter  increases  as  more  rarefied  drag  measurements 
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are  made  (e.  g. ,  see  Refs.  1,  4,  and  5).  If  the  experimental  scatter  is 
greater  than  the  effect  of  Tw,  or  if  sufficient  data  points  are  not  ob¬ 
tained  in  order  to  apply  a  statistical  analysis,  free-flight  data  are 
difficult  to  interpret  as  to  the  effect  of  Tw  only.  Furthermore,  the 
models  are  usually  small  and  made  of  lightweight  materials  having  low 
thermal  conductivity.  Thus,  the  accuracy  of  determining  their  mass 
and  the  effective  wall  temperature  experienced  by  the  molecules  de¬ 
creases.  In  view  of  the  advantages  and  disadvantages  of  the  balance 
and  free-flight  techniques,  and  because  a  suitable  drag  balance  was 
available  (Ref.  5),  the  balance  technique  was  chosen  for  this  investiga¬ 
tion. 


A  discussion  of  the  appropriate  dynamic  simulation  parameters  for 
hypervelocity  bluff  body  drag  coefficients  has  been  given  by  Potter  and 
Miller  (Ref.  6).  Because  some  of  the  data  to  be  considered  here  cor¬ 
respond  to  low  supersonic  Mach  numbers,  another  expression  for  cor¬ 
relating  sphere  drag  data  will  be  derived  and  its  validity  assessed  in 
conjunction  with  an  approximate  analysis  presented  in  Refs.  3  and  7. 
Comparisons  of  several  sets  of  experimental  data  with  the  approximate 
theory  will  be  made. 


SECTION  II 

EXPERIMENTAL  FACILITY  AND  NOZZLE  CALIBRATION 


The  10V  chamber  of  the  Aerospace  Division  of  VKF  is  a  cryo- 
genically  pumped  chamber  10  ft  in  diameter  and  20  ft  in  length.  It  has 
about  an  8-ft-diam  and  10-ft-long  internal  working  volume  when  the 
nominal  Mach  three  (M3)  or  Mach  six  (M6)  nozzle  is  installed  for  oper¬ 
ation  as  a  wind  tunnel.  The  principal  pumping  capacity  of  the  10V  for 
this  test  was  provided  by  620  ft 2  of  77°K  liquid -nitrogen  (LN2)  cryo- 
surfaces  and  240  ft ^  of  15  to  20°K  gaseous -helium  (GHe)  cryosurfaces. 
Although  8  kw  is  available,  only  a  4-kw  GHe  refrigeration  capacity  was 
used.  This  arrangement  permitted  the  continuous  pumping  of  8  to 
10  gm/sec  of  nitrogen  (or  argon)  at  a  chamber  pressure  of  about 
10' 4  mm  Hg. 

The  10V  instrumentation  and  hardware  arrangement  was  the  same 
for  this  test  as  described  in  Ref.  3.  A  schematic  of  the  chamber  is 
given  in  Fig.  1  (Appendix  I).  The  balance  used  to  measure  the  drag 
force  is  described  in  detail  in  Ref.  5.  Problems  associated  with 
balance  alignment,  vibration  damping,  thermal  control,  and  calibration 
(all  of  which  were  accomplished  remotely  in  the  10V  environment)  are 
also  discussed  in  Ref.  5. 
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■Drag  measurements  were  made  on  the  centerline,  10  in.-  down¬ 
stream  of  the  exit  plane  of  the  M3  nozzle.  The  flow  was  calibrated  at 
this  point  using  a  1-in.  -diam,  10-deg  internally  chamfered  pitot  probe. 

An'  MKS  Baratron®  (variable  capacitance  transducer)  was  used  to  measure 
both  the  total  pressure,  pQ,  and  the  pitot  pressure,  Pp.  Three  correc¬ 
tions  to  pitot  pressure  measurements  in  this  flow  regime  were  con¬ 
sidered:  (1)  thermal  transpiration  effect,  (2)  orifice  effect,  and 
(3)  viscous  correction.  The  thermal  transpiration  correction  was  found 
to  be  negligible.  Following  the  method  of  Kinslow  and  Potter  (Ref.  8) 
the  orifice  effect  was  found  to  contribute  less  than  one  percent  change 
in  measured  pressure  and  therefore  was  neglected.  The  external-flow 
viscous  correction  to  the  impact  pressure,  however,  was  a  large  as 
10  percent  and  therefore  not  negligible.  The  viscous  correction  data 
presented  in  Ref.  3  were  used  to  correct  the- pitot  measurements. 

These  data  were  curve  fitted  and  a  computer  program  written  to  carry 
out  the  necessary  iterations  for  data  reduction. 

The  .criterion  described  in  Ref.  3  for  merged  flow,  i.  e. ,  the 
boundary  layer  completely  filling  the  nozzle,  was  used  here.  The 
reservoir  Reynolds  number  based  on  throat  radius  (Re0)  r*.  see  Ref.  9) 
for  merged  flow  in  the  M3  nozzle  is  Re0j  r*  <  1000.  No  data  are  pre¬ 
sented  for  Re0j  r*  <  1000. 

The  axial  Mm  gradient  10  in.  downstream  of  the  conical  M3  nozzle 
is  about  0.  015  per  inch.  For  the  flow  conditions  of  these  measurements, 
this  corresponds  to  a  free-stream  dynamic  pressure  gradient  from  1.  2 
to  1.  3  percent  per  inch.  The  effect  of  this  gradient  on  for  a 
0.  25 -in.  -diam  sphere  was  therefore  negligible: 


SECTION  III  - 
DATA  ACQUISITION 


The  data  acquisition  procedure  was  somewhat  different  than  com¬ 
monly  used  to  make  drag  measurements  in  the  transition  regime. 
Usually  flow  conditions  are  established  and  force  measurements  are 
made  using  different  model  sizes,  thereby  obtaining  C;p  for  various  Kn,,, 
or  Re,,,  at  about  constant  Tw/T0O.  Or,  a  given  model  is  used  and  . 
pQ  varied  to  obtain  Cj)  versus  Kn^,.  Then,  other  flow  variables,  such 
as  TQ,  are  changed,  and  similar  data  are  acquired.  The  requirement 
of  large  Kn^,  or  small  Re^,  restricts  the  model  in  most  facilities  to  a 
size  such  that  Tw  cannot  be  varied  conveniently,  and  frequently  it  is 
difficult  to  measure  at  all.  Since  free-stream  mean  free  paths,  Xw,  on 
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the  order  of  an  inch  can  be  produced  in  the  10V,  sufficiently  large 
models  can  be  used  such  that  Tw  can  be  varied  and  continuously  moni¬ 
tored.  Therefore,  by  fixing  flow  conditions  and  model  size,  drag  data 
can  be  obtained  which  produce  vertical  lines  on  a  plot  of  Cj)  versus  Kn^, 
or  Re^,  since  each  data  point  is  associated  with  a  different  TW/T(D. 

The  standard  procedure  for  making  these  drag  measurements  was 
to  cool  the  sphere  to  some  temperature,  establish  flow  conditions,  and 
then  periodically  make  force  measurements  as  the  sphere  temperature 
increased  due  to  convective  heating  of  the  sphere  by  the  flowing  gas. 

Since  the  chamber  pressure  was  typically  about  10"?  mm  Hg  when  there 
was  no  flow  into  the  chamber,  there  was  insufficient  gas  density  to  de¬ 
crease  the  sphere  temperature  below  about  270°K  even  though  there  are 
77°K  and  15-20°K  cryopanels  in  the  chamber.  Also  radiation  heat  trans¬ 
fer  would  not  decrease  the  sphere  temperature  below  270°K  when  the 
sphere  was  located  in  the  center  of  the  chamber.  Therefore,  a  means  of 
convective  heat  transfer  was  supplied  to  the  chamber  by  a  1.  3  x  10"  ^ 
cc/sec  (standard)  calibrated  helium  leak  to  raise  the  chamber  pressure 
to  about  10" ^  mm  Hg.  Since  15°K  surfaces  will  not  effectively  pump 
helium,  the  helium  will  stay  in  the  chamber  until  removed  by  diffusion 
pumps,  and  thereby  the  medium  necessary  for  convective  cooling  of  the 
sphere  could  be  maintained  as  long  as  required.  Unfortunately,  because 
of  the  time  required  to  establish  flow  conditions  and  make  a  force  meas¬ 
urement,  the  lowest  Tw  at  which  data  were  taken  was  about  210°K.  The 
largest  value  of  Tw  attainable  was  limited  by  radiation  to  the  cryopanels. 
For  the  higher  T0  data  the  largest  Tw  attained  was  about  390°K.  Typical 
rates  of  increase  of  Tw  are  shown  in  Fig.  2. 

\ 

A  second  0.  25(±0.  0002) -in.  -diam  aluminum  sphere,  identical  to 
the  one  used  on  the  balance  sting  to  measure  the  drag,  was  used  to 
determine  the  sphere  wall  temperature.  This  second  sphere  was  sup¬ 
ported  by  the  leads  to  a  thermocouple  which  was  used  to  measure  the 
sphere  temperature.  It  was  attached  to  the  same  support  as  the  balance 
and  was  simultaneously  subjected  to  the  same  flow  conditions  as  the 
sphere  mounted  on  the  balance  sting.  From  the  time -dependent  solution 
for  heat  transfer  in  a  sphere  (Eef.  10),  the  temperature  at  the  thermo¬ 
couple  was  calculated  to  be  within  one  percent  of  the  surface  tempera¬ 
ture  in  less  than  two  seconds.  This  lag  was  therefore  negligible,  and 
the  thermocouple  reading  was  taken  as  Tw. 

The  method  used  to  reduce  and  determine  tare  forces  was  the  same 
as  that  described  and  illustrated  in  Ref.  3.  Tare  forces  for  these  meas¬ 
urements  varied  from  a  minimum  of  27  percent  for  the  T0  =  30G°K  data 
to  a  maximum  of  38  percent  for  the  T0  =  755°K  data.  An  example  of  the 
force  and  pitot  pressure  measurements  is  given  in  Fig.  3. 
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SECTION  IV 

EXPERIMENTAL  RESULTS 


Drag  measurements  were  made  at  three  different  reservoir  tem¬ 
peratures  of  nominally  300,  640,  and  750°K.  These  data  are  tabulated 
in  Table  I  (Appendix  II)  and  presented  in  the  form  of  C;q  versus  Re„  in 
Fig.(  4.  The  drag  data  definitely  illustrate  a  "grouping"  with  T0  when 
presented  in  this  form.  However,  it  should  be  noted  that  each  group 
of..T0  data  has  associated  with  it  a  different  range  of  flow  conditions, 
and  hence  each  T0  group  has  a  different  CDfm  which  it  is  assumed 

to  approach  asymptotically.  (CDfm  is  herein  based  on  complete  accom¬ 
modation  and  diffuse  reflection. )  For  these  data  the  M,,,  and  Tw/Tro 
ranges  associated  with  each  T0  group  maximize  the  difference  in  the 
CDfm  limit .  For  example,  the  T0  =  300°K  data  have  the  lowest  values 
of.M,,,  and  the  largest  values  of  Tw/Ta>,  each  contributing  to  increasing 
CDfm»  whereas  the  TQ  =  750°K  data  have  the  largest  values  of  M0  and 

lowest  values  of  Tw/Tm,  each  contributing  to  decreasing  CDfm-  There¬ 
fore,  it  is  appropriate  to  consider  the  normalization  of  the  drag  data  by 
the  CDfm  value  associated  with  each  data  point.  The  results  of  such  a 
normalization  will  be  presented  subsequently. 

Obviously,  the  normalization  of  Cd  by  CDfm  is  sufficient  for  cor¬ 
relating  free -molecular  drag  data  if  the  theory  used  to  calculate  CDfm 
is  valid.  However,  these  drag  data  correspond  to  Kn^  of.  order  one  and 
therefore  are  subject  to  the  physics  of  the  more  complicated  transition 
regime.  Since  no  general  solution  has  been  found  for  the  Boitzmann 
equation  which  describes  the  flow  about  a  body  in  this  regime,  the  appro¬ 
priate  parameters  for  correlating  these  drag-data  are  not  available.. 
Fortunately,  approximate  solutions  of  the  Boltzmann  equation  and  solu¬ 
tions- to  approximate  models  of  transition  flow  have  exposed  useful 
parameters.  A  summary  of  most  of  the  parameters  applicable  to  hyper¬ 
sonic  flows  is  given  by  Potter  and  Miller  (Ref.  6).  Attention  will  be 
given  here  to  correlating  supersonic  drag  data.  A  new  correlation 
parameter  will  be  derived  and  its  validity  investigated  by  using  it  to 
correlate  these  as  well  as  other  drag  data.  Although  the  purpose  of  this 
report  is  not  to  derive' and  assess  various  correlation  parameters,  the 
use  of  a  parameter  to  correlate  these  data  is  important  since  variables 
which  can  affect  sphere  drag  in  this  regime  change  for  each  data  point 
in  the  present  experiment. 
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SECTION  V 

CORRELATION  PARAMETER 


The  Knudsen  number,  Knw  =  Aw/d,  based  on  the  mean  distance 
traveled  by  molecules  emitted  from  the  surface  of  a  body  before  experi¬ 
encing  a  collision  with  free -stream  molecules,  has  been  shown  in 
various  theoretical  sphere  drag  investigations  (Refs.  7  and  11  through  16} 
to  be  of  importance.  The  derivations  of  Aw  have  assumed  that  the 
emitted  molecules  travel  at  some  average  velocity,  and  in  some  in¬ 
stances  Aw  is  restricted  to  cold-wall,  high-Mach-number  conditions.  It 
is  of  interest,  therefore,  to  consider  fundamentally  the  derivation  of  Aw 
and  to  determine  whether  or  not  a  more  accurate  expression  for  Aw  can 
be  obtained  without  formidable  mathematical  difficulties. 


Consider  two  classes  of  smooth,  rigid,  elastic,  spherical  molecules 
denoted  as  1  and  2.  For  such  molecules  no  problem  exists  concerning 
the  concept  of  a  free  path  since  the  molecules  affect  each  other's  motion 
only  at  collisions.  Following  the  development  of  Chapman  and  Cowling 
(Ref.  17),  the  total  number  of  collisions  occurring  per  unit  volume  and 
unit  time  between  pairs  of  molecules  of  class  1  with  class  2  is  given  by 


where 


and 


N12  -////;/  rrg  °12  ^1  ^2  duj  dvj  dwj  du2  dv2  dw2 


g  =  C(u x  -  U2)2  +  (vj  -  V2)2  +  (wj  -  w2)2]^ 


(712  =  —  +  a2) 


(1) 

(2) 

(3) 


In  Eq.  (3),  and  are  the  diameters  of  the  molecules  of  class  1  and 
class  2,  and  f^  and  f2  are  the  distribution  functions  of  class  1  and 
class  2.  For  the  coordinate  system  shown  in  Fig.  5,  becomes, 
after  nondimensionalizing, 


2jt  C712  Iij  n2  V  2RTi 


h  x 


e-[(6-Si  sin 0)2  +  O?l  +s,  COS^2+^2j-(^2+722  +  C22)d^idj?id<id&dj?2d<2  (4) 


where 
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yj~2 RT7 


li  =  Kfl— r^)2  +  (t?1  -n/2)2  +  (^1  -  r^2)2^  2  .  (5e) 

r  =  (J2/Ti)V3  <5f> 

and.i  may  take  on  the  values  1  or  2.  Note  the  additional  2  as  a  coeffi¬ 
cient  in  the  integrand  of  Eq.  (4).  This  2  is  required  in  order  to  satisfy 
continuity  since  the  limits  of  §2  are  "®  to  0,  and  hence  only  one  half  of 
the  molecules  emitted  from  the  surface  element  in  Fig.  5  would  other¬ 
wise  be  accounted  for.  The  mean  distance  traveled  by  molecules  of 
class  2  between  successive  collisions  with  those  of  class  1  is  given  by 
(cf.  Ref.  17) 

.  1*2  v2 

X21  =  - 


where  v 2  is  the  mean  speed  of  class  2  molecules.  Taking 
v2'=  (9ttRT2/8)1/2  (Ref.  13),  Xn  =  (>f2_7rn1a112)'1  (Ref.  17),  and 
ai  =  ?2  <since  the  same  size  molecules  are  being  considered  and  no 
dependence  of  u  on  temperature  is  assumed  at  this  point),  one  obtains 

*21  _  Syfin  rr3  r 
An  8  I 

where 

1=  /  f  /  /  f  f  h  sin/3)Z+(j;i +Sicos^)2 -!-Ci2  +  ^22 +»/22  +  C22] 

--  -V  -«  i  Jo  „  ^  ,  _. 

'  dfidr;  id^i  ds2tl7?2^C2  \8) 

and  h  is  given  by  Eq.  (5e). 


The  integral  I  given  by  Eq.  (8)  has  not  been  integrated  in  closed 
form.  Approximate  results  have  been  obtained  for  Eq.  (7)  by  replacing  g 
in  Eq.  (1)  by  the  relative  mean  speeds  of  the  class  1  and  2  molecules, 
see,  e.  g. ,  Ref.  18.  The  resulting  expression  has  been  used  for  the 
mean  free  path  of  emitted  molecules  from  a  surface  in  several  theo¬ 
retical  investigations,  e.  g. ,  Refs.  7,  11,  13,  and  15.  Of  course  the 
accuracy  of  this  approximation  can  only  be  determined  by  comparison 
with  exact  solutions  of  Eq.  (8).  Even  numerical  solutions  to  Eq.  (8)  are 
rather  difficult  and  would  require  a  good  deal  of  computer  time  for 
accurate  results  for  the  several  values  of  t,  and  3  which  should  be 
investigated.  In  view  of  this,  it  is  of  interest  to  present  anotner  approxi¬ 
mate  solution  to  Eq.  (8)  which  will  be  shown  to  be  in  reasonable  agreement 
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with  the  numerical  solutions  of  Eq.  (8)  which  have  been  obtained.  Also, 
this  approximate  solution  will  provide  an  upper  bound  for  ^-21  An 
(Eq.  (7)),  whereas  the  approximate  result  given  by  using  mean  speeds 
for  g  in  Eq.  (1)  provides  a  lower  bound  for  ^-21  All  when  compared  to 
the  numerical  solutions. 

Consider  the  integrand  of  Eq.  (8).  Notice  that  it  is  non-negative 
and  attains  its  larger  values  at  relatively  small  values  of  independent 
variables  ?2»  ^2*  anc*  ^2  because  of  the  exponential  term.  Because 

of  the  terms  sin  0)  and  (r?i  4  cos  0)  in  the  exponential,  the 

integrand  will  attain  its  larger  values  for  near  (S^  sin  0)  and  near 
(-S^  cos  0).  It  appears  that  in  rarefied  flow,  molecules  emitted  near  0 
of  it  12  have  more  effect  on  sphere  drag  than  those  emitted  at  locations 
remote  from  the  stagnation  point.  Therefore,  Eq.  (8)  will  be  considered 
‘  for  the  case  |  Si  sin  /3 1  »  |  —  S |_  cos  0 1 .  With  | S i  sin  /3  |  »  |  -  Si  cos  0| 
the  integrand  of  Eq.  (8)  is  larger  when  ?i  is  on  the  order  of  Si  sin  0, 
because  ?i  contributes  more  to  increasing  h  before  the  exponential  term 
dominates  than  do  the  other  independent  variables.  The  first  term  in  h 
(Eq.  (5e))  is  therefore  retained,  and  h  is  approximated  by  h  =  £i  -  %27' 
Of  course  £2  will  be  small  compared  to  ?i,  when  the  integrand  of  Eq.  (8) 
is  very  large  (because  the  exponential  term  dominates  when  ?2  is  large) 
and  hence  will  contribute  about  the  same  as  the  terms  dropped.  How¬ 
ever,  notice  that  the  integral  of  (-?2T)  times  the  exponential  is  greater 
than  zero  because  only  negative  values  of  £2  are  considered.  Therefore, 
since  the  contribution  of  the  integral  of  the  terms  which  were  dropped  is 
greater  than  zero,  then  the  approximation  is  slightly  more  accurate  if 
(-£2t)  is  retained.  The  resulting  integral  is 

'-rrrrrxv*** 

e-U6-Si  sin0)2  +  (7?1tS1  cos/3)2  +  Ci2  +  ^2+V22j-C22}^li7]li(rl^2ir]2^2  (9) 


There  is  no  problem  in  integrating  Eq.  (9),  see,  e.  g. ,  the  Appendix 
of  Ref.  19.  The  result  is 


I  =  [[1+  erf  (Si  sin/3)]  [\0TSisin0  +  r]  +  e.-Sl2  sin20j 


(10) 


Substituting  Eq.  (10)  into  Eq.  (7)  for  I  gives 


^21 

*11 


3  n't! 


2{[H-erf  (Sx  sid/3)]  [/ V  v  Si  sin/3  +  r]  +  e"Sl2  S'™ZP] 


(ID 


8 


AEDC-TR-71-83 


Therefore,  for  Si  sin  /3  of  order  one  and  larger,  Eq.  (11)  can  be  rather 
accurately  approximated  by 

A21  _  3  n\ 2  r _ 

Xu  4{>|'F‘  Si  sin/3  +  t)  (12) 

This  result  is  applied  to  the  problem  at  hand  by  considering  the 
class  1  molecules  to  be  free-stream  molecules,  =  Am,  and  the 
class  2  molecules  to  be  molecules  emitted  from  the  surface  of  the 

sphere,  A21  =  Aw.  Then  t  becomes  V  Tw/Tm.  Using  Sj  =  S^, 

Sw  =  S^/  \J  TW/T(D  ,  and  the  approximation  sin  /3  =  1,  one  obtains 

Aw  _  3  n\fi 

Am  4<V5FSw  +  1)  (13) 

The  adequacy  of  Eq.  (13)  is  now  assessed  by  comparisons  with 
results  using  numerical  solutions  for  I  in  Eq.  (7).  In  view  of  Eq.  (13), 
it  .appears  more  convenient  to  consider  the  reciprocal  of  Aw/A0o  since 
at  least  the  approximate  solution  is  linear  in  Sw.  Results  obtained 
numerically,  by  the  present  approximation,  and  by  using  mean  speeds 
for  g  in  Eq.  (1)  are  all  presented  in  Fig,  6.  The  numerical  solutions 
in  Fig.  6  were  obtained  for  1/2  <  t  <  5  which  corresponds  to 
1/4  <  Tw/T,,  <  25  for  each  Sw. 

As  mentioned  previously,  the  approximate  solutions  provide  upper 
and  lower  bounds  to  the  numerical  solutions  for  these  values  of  the 
parameters  t  and  S,,,.  It  is  of  interest  to  note  that  the  slopes  of  the 
approximate  results  are  identical  and  correspond  to  about  the  mean 
slope  of  the  numerical  solutions  for  S,,,  >  2.  Equation  (13)  is  in  better 
agreement  with  the  numerical  results  as  becomes  large,  and  since 
this  regime  is  of  more  interest  here  (and  in  the  approximate  theo¬ 
retical  analysis  where  the  other  result  has  been  used),  it  seems  more 
appropriate  to  use  Eq.  (13),  i.  e.,  take 

Aw  3ir\/2  Kn«j  Kno- 

r  2<vF  Sw  +  1)  \fn  Sw  +  1  (  14) 


It  should  be  noted  that  the  result  given  in  Eq.  (13)  is  related  to  that 
of  Ref.  13,  i.  e. , 


( Aw/Aoo)  Ref<  l3 

(  Aw/A«,)Eq>  (l3) 


Sw  y!7  +  1 
S +  2.35  • 


1  when  Sw\/F »  1 


Therefore,  even  though  either  approximation  may  be  used,  it  has  been 
found  that  use  of  Eq.  (13)  results  in  slightly  better  correlation  of 
experimental  data. 
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SECTION  VI 

THEORETICAL  AND  EXPERIMENTAL  COMPARISONS 


The  more  appropriate  available  sphere  drag  data  with  which  to  com¬ 
pare  the  present  results  are  those  of  Potter  and  Miller  (Ref.  6),  Davis 
and  Sims  (Ref.  20),  Smolderen,  et  al.  (Ref.  21),  and  Whitfield  and 
Stephenson  (Ref.  3).  Also,  the  lower  density  data  of  Bailey  and  Hiatt 
(Ref.  1)  are  considered.  The  data  from  these  sources  are  compared 
with  the  present  results  in  Figs.  7  through  10  in  terms  of  several  param¬ 
eters.  Figures  7  and  8  present  the  data  in  terms  of  the  free-stream 
Reynolds  and  Knudsen  numbers.  Owing  to  the  variations  in  Mm  and 
Tw/Tas,  correlation  of  the  results  would  not  be  expected.  In  Fig.  9 
these  data  are  presented  with  Cq  normalized  by  CDfm»  where  CDfm  is 

calculated  by  using  the  particular  S,,,  and  Tw/T w  values  associated  with 
each  data  point.  A  further  improvement  is  illustrated  in  Fig.  10  by 
using  Kn^,  /( yflT  Sw  4  1),  from  Eq.  (14),  in  place  of  Kn^.  The  approxi¬ 
mate  theory  in  Fig.  10  was  obtained  by  using  Xw/r  given  by  Eq.  (14)  in 
the  expression  derived  for  the  sphere  drag  coefficient  in  Ref.  7.  The 
resulting  equation  is 


where 


Cp 

CDf, 


l 


K  ,  _  Kn, 


\rr  Sw  +  1 


(15) 


(16) 


The  comparison  between  experiment  and  theory  and  experiment  in 
Fig.  10  is  reasonable  for  >  10 _1  considering  that  the  experimental 

results  consist  of  data  over  an  M,,,  range  from  2.  5  to  8.  3  and  a  Tw/T0 
range  from  0.  8  to  15.  Also,  the  measurements  made  by  Davis  and  Sims 
(Ref.  20)  were  taken  in  a  carbon  dioxide  (CO2)  gas,  whereas  the  other 
measurements  were  made  in  nitrogen  and  air.  Although  the  range  data 
of-  Bailey  and  Hiatt  (Ref.  1)  do  not  overlap  the  conditions  of  the  present 
measurements,  the  comparison  in  Fig.  10  illustrates  that  their  data 
are  in  agreement  with  the  data  of  Potter  and  Miller  (Ref.  6),  Davis  and 


Sims  (Ref.  20),  and  the  theory;  all  three  of  which  are  in  good  agreement 
and  represent  about  the  mean  drag  value  of  the  present  data. 
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Since  slightly  more  accurate  numerical  solutions  were  obtained  for 
^■w/^oo  (Fig.  6),  it  is  of  interest  to  determine  the  changes  attributable 
to  using  numerical  values  of  Xy^/X^  in  place  of  Eq.  (13).  In  order  to  do 
this,  a  mean  linear  curve  was  placed  through  the  appropriate  range  of 
conditions  in  Fig.  6  corresponding  to  the  data  in  Fig.  10,  and  the  equa¬ 
tion  of  the  linear  curve  was  determined.  As  a  result,  X^/X^  was  pro¬ 
portional  to  (Sw  4  1)  and  therefore  Aw/r  -  Kn^/fSw  4  1).  Figure  11 
represents  the  correlation  provided  by  the  numerical  solutions  for  the 
same  data  considered  in  Figs.  7  through  10.  The  approximate  theory 
was  obtained  by  substituting  the  linear  equation  for  Aw/r  obtained  from 
the  numerical  solutions  into  the  sphere  drag  expression  of  Ref.  7.  Close 
examination  of  Fig.  11  indicates  that  little  change  in  the  correlation  was 
obtained. 


A  more  critical  comparison  between  theory  and  the  variation  of  the 
present  drag  data  with  Tw  is  provided  by  Fig.  12.  Data  used  for  this 
comparison  are  those  taken  at  essentially  constant  and  Kn^,  and 
therefore  any  change  in  drag  should  be  due  only  to  a  change  in  Tw.  For 
constant  and  Kn^,  the  theoretical  Cj)  of  Ref.  7  depends  only  on 
Tw/TOJ  and  hence  curves  of  constant  and  Kn^  can  be  constructed 
for  comparison  in  Fig.  12.  The  experimental  data  and  theory  have 
almost  identical  slopes  in  Fig.  12,  and  the  comparison  is  considered 
good.  Also,  the  variation  of  the  present  data  and  theory  with  T^/T^  in 
Fig.  12  is  further  substantiated  by  the  data  of  Potter  and  Miller  (Ref.  6), 
Davis  and  Sims  (Ref.  20),  and  Bailey  and  Hiatt  (Ref.  1).  The  variation 
of  the  present  drag  data  with  Sw  is  compared  with  theory  in  Fig.  13. 


Calculated  values  of  Cq  from  Eq.  (15)  (at  constant  with  7  =  1.  4 
and  Tyf/T^  =  1)  are  presented  as  a  function  of  Re,,,  in  Fig.  14.  With 
7  =  1.4  and  T-w/T^  =  1,  Eq.  (15)  can  be  written  in  terms  of  M0  and  Re^ 
for  M.  >  2.5  and  =  (7T/2RTJ1  >2  palpn.  viz.. 


2.0004 


CD  = 


1.412 

Moo 


2,857 

Mm2 


1  4 


('  Moo  \2 

,  +  _ _^L) 

L  1.483  Mcc  4-  1/ 


(17) 


This  Cd  is  not  valid  for  continuum  flow  since  the  drag  becomes  in¬ 
sensitive  to  Tw  at  high  densities,  and  CD/CDfm  are  usually 

something  other  than  one -half  as  Kn^  -►  0.  However,  it  has  been 
shown  to  be  valid  through  the  transition  regime  to  the  slip  flow  regime 
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(Re^  22).  In  terms  of  limits,  Eq.  (15),  and  hence  Eq.  (17),  can  be  used 
for  Kw>  d  =KnB/(v^ Sw  +  1)  >  0.  1  or  Kr<  ^  =  Kn0/(^  Sr  +  1)  >  0.  1.  For 
the  conditions  of  the  present  experiment  atM,,*  3,  this  implies  Re^  <  0(10). 

1  Although  hypersonic  cold-wall  data  are  not  of  prime  importance  here, 
they  become  of  interest  when  considered  for  correlation  of  data  such  as 
used  in  Figs.  7  through  10.  A  plot  of  several  sets  of  high  enthalpy  (HQ) 
hypersonic  data  (Refs.  6,  13,  21,  and  23  through  25)  in  terms  of  KWj  d 
(Fig.  15,  open  symbols)  reveals  that  the  cold- wall  data  have  larger 
values  of  Co/CDfm.  at  same  values  of  Kw>  d  than  the  lower  and/or 

lower  H0  data  considered  in  Figs.  7  through  10.  Although.it  is  only 
speculation  at  this  point,  an  explanation  might  be  that  the  energy  of  the 
gas  molecules  on  the  outermost  surface  of  the  sphere  is  rapidly  in¬ 
creased  above  that  corresponding  to  Tw  by  collisions  with  the  high  energy 
free-stream  molecules.  If  this  is  the  case,  an  increase  in  drag  occurs 
because  of  the  higher  momentum  (pressure)  exerted  on  the  body  by  the 
molecules  which  leave  the  surface.  The  work  of  Kinslow  and  Potter 
(Ref.  8)  does  not  address  this  problem  specifically;  however,  it  is  con¬ 
sidered  in  search  of  an  estimate  of  this  effect. 


Kinslow  and  Potter  (Ref.  8)  used  a  two-sided  Maxwellian  distribu¬ 
tion  function  to  describe  the  gas  in  the  vicinity  of  a  surface.  For  a  bluff 
body  stagnation  region  their  expression  for  the  surface  pressure  pw  is 
approximated  by 


P  w 

Pifm 


(18) 


where  Pifm  =  mn'RT"(Tw/T')^^  is  the  pressure  in  a  cavity  in  the  sur¬ 
face  where  the  cavity  orifice  is  much  smaller  than  the  mean  free  path  in 
the  cavity,  and  the  reference  temperature  Tr  is  defined  as 


Tr  = 


T0  +  T+ 
2 


(19) 


The  minus  superscript  denotes  incoming  molecules  and  the  plus  super¬ 
script  denotes  outgoing  molecules.  The  present  approximation  made  in 
Eq.  (18)  is  the  use  of  the  exponent  1/2  instead  of  the  more  accurate  value 
which  may  vary  between  approximately  1/2  and  1/(2.  3).  Since  Pifm  is 

proportional  to  (Tw)i^,  the  actual  surface  pressure  pw  is  given  approxi¬ 
mately  by  replacing  Tw  with  Tr  in  Pifm- 


It  seems  appropriate  to  investigate  the  change  in  the  normalized 
drag  data  caused  by  using  Tr  in  place  of  Tw  in  determining  the  pressure 
attributable  to  emitted  molecules  in  calculating  CDfm.  The  reference 


12 


AEDC -TR-71-83 


temperature,  Tr,  was  used  to  calculate  both  QDfm  and  Kr>  ^  (i.  e. , 

Rr<  d  =  Kn a,/(^V" Sr  4  1))  in  Fig.  15  (closed  symbols).  The  solid  curve 
in  Fig.  15  is  given  by  Eq.  (15).  As  a  result,  for  the  data  considered 
in  Figs.  7  through  10,  Tr  =  Tw,  and  quite  good  correlation  of  all  the 
data  is  obtained  by  using  Tr.  Moreover,  the  scatter  in  the  data  of  a 
given  set  is  decreased  somewhat;  note  particularly  the  17.  43  <  <  21.00 

data  of  Ref.  6. 

Obviously,  the  argument  presented  for  the  use  of  Tr  for  correlating 
drag  data  is  not  rigorous.  At  this  point  the  only  justification  for  its  use 
is  that  apparently  it  works.  It  should  be  pointed  out  that  Kinslow  and 
Potter  (Ref.  8)  showed  that  the  existing  data  associated  with  the  so- 
called  orifice  effect  can  be  correlated  using  Tr  to  calculate  a  Knuds en 
number  based  on  the  orifice  diameter.  The  physics  behind  the  usefulness 
of  Tr  for  correlating  drag  data  should  be  investigated  further. 


SECTION  VII 
CONCLUSIONS 


•  1  The  present  drag  data  for  constant  T0  are  correlated  quite  well  by 
the  parameter  used  in  Fig.  10  (the  symbol  and  T0  correspondence  is 
identified  in  Fig.  4);  however,  there  is  a  two-  to  four-percent  differ¬ 
ence  in  the  drag  coefficients  of  the  three  groups  of  constant  T0  data. 
Probably  the  most  plausible  explanation  for  this  discrepancy  is  the 
presence  of  small  errors.  Possible  errors  which  are  considered  to 
contribute  less  than  one-percent  change  in  the  results  are  (1)  thermal 
transpiration  correction  to  measured  pQ,  (2)  error  involved  in  meas¬ 
uring  T0,  (3)  zero  shift  of  drag  balance  caused  by  thermal  environment, 
and  (4)  inaccuracy  in  positioning  the  balance  and  pitot  probes.  A  more 
likely  source  of  error  is  introduced  by  the  actual  flow  calibration.  The 
necessary  corrections  to  the  pitot  pressure  measurements  increase 
with  T0  to  a  maximum  of  about  10-percent  correction'at  T0  =  755°K. 

It  is  entirely  possible  that  a  combination  of  the  smaller  errors  men¬ 
tioned  above,  errors  in  actually  measuring  pitot  pressures,  and  errors 
introduced  by  the  corrections  applied  to  the  pitot  pressures  could  account 
for  the  two-  to  four-percent  variation  in  among  the  three  T0 

groups. 

It  should  be  kept  in  mind,  however,  that  the  total  deviation  of  the 
present  data  is  on  the  order  of  the  scatter  in  existing  free -flight  data  in 
this  regime  acquired  supposedly  at  constant  flow  conditions.  Moreover, 
the  mean  of  all  the  present  data  is  shown  in  Fig.  10  to  be  in  good 
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agreement  with  available  experimental  data.  The  scatter  among  the 
presented  data  at  constant  T0  is  less  than  ±1  percent  (Ref.  3),  and  the 
variation  or  slope  of  Cd  with  Tw  (Fig.  12)  is  considered  accurate  with¬ 
in  one  percent. 

The  parameter  provided  by  Eq.  (14)  is  useful  in  correlating  a  class 
of  sphere  drag  data  for  Kw<  ^  >  0.  01  and  low  total  enthalpy.  Equation  (15) 
is  shown  in  Fig.  10  to  be  in  rather  good  agreement  with  the  data  of  prime 
interest. 
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At,  minutes 

Fig.  2  Change  in  Sphere  Wall  Temperature  with  Time 


20 


AEDC-TR-71-83 


0  0.1  0.2  0.3 

PQ,  mm  Hg 

Fig.  3  Force  and  Pitot  Pressure  Measurements  10  in.  Downstream  of  Mach  3  Nozzle  Exit 


Fig.  4  Summary  of  Experimen 
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Fig.  5  Coordinate  System 
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Fig.  6  Numerical  and  Approximate  Analytical  Solutions  for  \J'KW 
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Fig.  8  Sphere 
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Fig.  10  Sphere  Drag  Data  in  Terms  of  CD/CD,  and  Kn 
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Fig.  11  Sphere  Drag  Data  According  to  Correlation  Provided  by  Numerical  Solutions  for  Xw/X. 
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Fig.  12  Effect  of  Wall  Temperature  on  Sphere  Drag 
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Fig.  13  Variation  of  Present  Drag  Results  with  Wall  Speed  Ratio 
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TABLE  I 

SPHERE  DRAG  AND  FLOW  CONDITIONS  IN  THE  MACH  3  NOZZLE 


Re* 

Kn* 

Tw/T„ 

sw 

' 

CD 

cD/cDfm 

p0,  mm  Hg 

T 

^0* 

°K 

Tw,  °K 

R-w,  d 

3,  450 

9.  274 

0. 5546 

2.  753 

1.  740 

2.406 

0. 8263 

0. 2000 

300.  0 

244. 27 

0.  1358 

3.  441 

8.  947 

0. 5734 

2.  368 

1.  871 

2.  327 

0. 8119 

0.  1920 

210.94 

0.  1329 

3.  445 

9.  070 

0. 5662 

2.  559 

1.  802 

2.  343 

0.  8108 

0.  1950 

227.  60 

0. 1350 

3.445 

9.  070 

0.  5662 

2.  684 

1.  759 

2.  360 

0.  8123 

0. 1950 

238. 72 

0. 1375 

3.445 

9.  070 

0. 5662 

2.  778 

1.  729 

2.  388 

0.  8186 

0. 1950 

247.  05 

0. 1393 

3.446 

9.  Ill 

0. 5638 

2.841 

1.  710 

2.  417 

0. 8265 

0.  1960 

252. 60 

0. 1399 

3.  383 

7.  167 

0. 7038 

2.  769 

1.  701 

2.  491 

0.  8484 

0. 1490 

252.  60 

0.  1753 

3.  338 

6.  112 

0.  8143 

2.  725 

1.  692 

2.  498 

0.  8478 

0. 1240 

253. 16 

0. 2036 

3.  308 

5.  514 

0. 8944 

2.  655 

1.  698 

2.557 

0.  8671 

0. 1100 

249. 83 

0.  2230 

3.  319 

5.  728 

0. 8639 

2.  668 

1. 700 

2.517 

0.  8545 

0. 1150 

249. 83 

0.  2152 

3.  330 

5.  942 

0. 8356 

2.  674 

1.  704 

2.  534 

0.  8610 

0.  1200 

249. 27 

0. 2078 

3.  352 

6.409 

0. 7797 

2.  698 

1.  707 

2.511 

0.8547 

0. 1310 

249. 27 

0.  1937 

3.  372 

6.  873 

0. 7314 

2.  726 

1.  709 

2.  526 

0.8605 

0. 1420 

249. 83 

0.  1816 

3.  385 

7.  208 

0. 7000 

2.  747 

1.  709 

2.487 

0.  8480 

0.  1500 

250. 38 

0. 1738 

3.400 

7.  626 

0. 6647 

2.  764 

1. 711 

2.509 

0.  8563 

0. 1600 

250. 38 

0. 1648 

3.414 

8.  040 

0.  6330 

2.  792 

1.  709 

2.  470 

0.  8433 

0. 1700 

251. 49 

0.  1571 

3.428 

8.  495 

0. 6016 

2.821 

1. 708 

2.449 

0.  8364 

0. 1810 

252. 60 

0.  1494 

3.439 

8.  865 

0. 5783 

2.  834 

1. 709 

2.  439 

0.  8338 

0. 1900 

252. 60 

0. 1435 

3.450 

9.  274 

0. 5546 

2.847 

1.  711 

2.  430 

0.  8311 

0. 2000 

252. 60 

0. 1375 

3.  593 

6.  571 

0.  8152 

1.  683 

2.  317 

2.  188 

0.8029 

0.  4400 

641.  5 

301. 49 

0. 1596 

3.594 

6.  627 

0.  8085 

1.  768 

2.  261 

2/207 

0.8062 

0.4440 

641.5 

316. 49 

0. 1614 

3.594 

6.  655 

0.  8053 

1.  877 

2.  195 

2.  245 

0.  8153 

0. 4460 

641.  5 

335. 94 

0. 1646 

3.  586 

6.  157 

0. 8685 

1.964 

2.  141 

2.  316 

0.  8367 

0. 4130 

644.  3 

354. 27 

0. 1811 

3.587 

6.  171 

0. 8666 

2.  057 

2.  092 

2.  360 

0.  8486 

0.  4140 

644.  3 

370.  94 

0. 1840 

3.586 

6.  157 

0. 8685 

2.  143 

2.  050 

2.  389 

0.8557 

0. 4130 

644.  3 

386. 49 

0. 1875 

3.550 

4.  474 

1.  183 

1.  417 

2.  495 

2.  264 

0.  8403 

0. 2870 

633.  2 

254. 83 

0.2182 

3.551 

4.  532 

1.  168 

1.  668 

2.  301 

2.  397 

0.  8769 

0. 2910 

633.  2 

299. 83 

0. 2301 

3.552 

4.  562 

1.  161 

1.  730 

2.  259 

2.  406 

0.  8770 

0.  2930 

633.  2 

.  310.94 

0.  2320 

3.552 

4.  562 

1.  161 

1.  823 

2.  201 

2.  452 

0.8895 

0. 2930 

633.  2 

327.  60 

0. 2368 

3.623 

5.  034 

1.073 

1.  327 

2.  632 

2.  077 

0.  7807 

0. 4240 

755.  6 

276. 49 

0. 1894 
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TABLE  I  (Concluded) 


He* 

Kn,,, 

H 

h" 

e- 

sw 

cD 

<:D/CDtm 

pQJ  nim  II g 

T0,  °K 

Tw,  °K 

Kw>  d 

3.  623 

5.  023 

1.075 

1.  438 

2.  527 

2.  191 

0.  8179 

0. 4230 

755.  6 

299. 83 

0.  1962 

3.  626 

5.  177 

1.  044 

1.  526 

2.456 

2.  146 

0.  7970 

0.  4370 

317.  60 

0.  1951 

3.  627 

5.  209 

1.  038 

1.  612 

2.  390 

2.  181 

0.  8060 

0.  4400 

335.  38 

0.  1982 

3.  628 

5.  242 

1.  032 

1.  663 

2.  354 

2.  186 

0. 8057 

0.  4430 

345. 94 

0.  1995 

3.  613 

4.  622 

1.  166 

1.  268 

2.  684 

2.  157 

0.  8132 

0. 3870 

265. 38 

0. 2024 

3.  615 

4.  711 

1.  144 

1.  373 

2.  581 

2.  174 

0.  8140 

0. 3950 

287. 05 

0. 2052 

3.  616 

4.  734 

1.  139 

1.  442 

2.  519 

2.  196 

0.  8190 

0.  3970 

301.  49 

0. 2084 

3.  616 

4.  756 

1.  134 

1.528 

2.  448 

2.  225 

0.  8254 

0. 3990 

319.  27 

0. 2123 

3.  617 

4.  779 

1.  128 

1.  613 

2.  382 

2.  238 

0.  8264 

0.  4010 

337.  05 

0. 2160 

3.  617 

4.  779 

1.  128 

1.  640 

2.  363 

2.  251 

0.  8298 

0.  4010 

342.  60 

0. 2175 

3.  617 

4.  790 

1.  126 

1.  714 

2.  31 1 

2.  266 

0.  8321 

0.  4020 

358.  16 

0. 2209 

3.  617 

4.  790 

1.  126 

1.  778 

2.  269 

2.  284 

0.  8357 

0. 4020 

371. 49 

0. 2242 

3.  603 

4.  165 

1.  200 

1.401 

2.  547 

2.  235 

0.  8346 

0.  3450 

752.  8 

293.  16 

0. 2339 

3.  604 

4.  188 

1.  283 

1.451 

2.  503 

2.  257 

0.  8403 

0.  3470 

303.  72 

0. 2360 

3.  604 

4.  211 

1.  276 

1.  547 

2.  424 

2.  283 

0. 8453 

0.  3490 

323.  72 

0. 2409 

3.  603 

4.  176 

1.  286 

1.  698 

2.  314 

2.  315 

0.  8497 

0. 3460 

355. 38 

0. 2522 

3.  617 

4.  836 

1.  1 15 

1.  227 

2.  732 

2.  061 

0.  7794 

0. 4040 

255.  38 

0. 1908 

3.  618 

4.  892 

1.  103 

1.  334 

2.  621 

2.  115 

0.  7943 

0. 4090 

277.  60 

0.  1953 

3.  619 

4.  961 

1.  087 

1.  452 

2.513 

2.  151 

0.  8018 

0. 4130 

750.  0 

300.  94 

0.  1994 

3.  620 

4.  995 

1.  080 

1.  528 

2.  450 

2.  163 

0.  8028 

0.  4160 

316.49 

0. 2022 

3.  620 

5.  006 

1.  078 

1.  622 

2.  378 

2.  191 

0. 8089 

0.  4170 

335.  94 

0. 2067 

3.  620 

5.  017 

1.  076 

1.  724 

2.  307 

2.  216 

0.  8134 

0. 4180 

357.  05 

0.  21 14 

3.  620 

5.  017 

1.  076 

1.  759 

2.  284 

2.  232 

0.  8177 

0. 4180 

364. 27 

0.  2131 

3.  624 

5.  125 

1.  054 

1.  204 

2.  764 

2.  051 

0.  7771 

0.4300 

249. 83 

0.  1787 

3.  626 

5.  191 

1.  041 

1.  253 

2.  710 

2.  067 

0.  7807 

0.  4360 

752.  8 

259. 83 

0.  1794 

3.  627 

5.  257 

1.  029 

1.  401 

2.  564 

2.  120 

0.  7935 

0. 4420 

290.  38 

0.  1855 

3.  620 

5.  301 

1.  021 

1.  501 

2.  478 

2.  147 

0.  7987 

0. 4460 

310. 94 

0.  1893 

3.  629 

5.  312 

1.  019 

1.  568 

2.  42b 

2.  164 

0.  8020 

0. 4470 

324.  83 

0.  1923 

3.  629 

5.  232 

1.  017 

1. 635 

2.  374 

2.  170 

0.  8011 

0.4480 

338.  72 

0. 1952 

3.  633 

I _ 

5.  426 

0.  0982 

1.  316 

2  650 

2.  072 

0.  7800 

0.  4600 

755.  6 

273. 16 

0. 1752 
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TABLE  II 

SPHERE  DRAG  AND  FLOW  CONDITIONS  IN  THE  MACH  6  NOZZLE 


Re0 

Knm 

Tw/T. 

S 

w 

CD 

ca/c%m 

PQ.  mm  Hg 

T 

‘o* 

•K 

Tw  °K 

Kw,  d 

7.  223 

2.  944 

3.  658 

4.  087 

2.  989 

2.  446 

0. 9984 

1.514 

80 

[).  0 

286.  00 

0.  5808 

7.  383 

3.  656 

3.  011 

4.  395 

2.  947 

2.  432 

0.  9913 

1.  991 

295. 50 

0.  4839 

7.  504 

4.  432 

2.  524 

4.  689 

2.  899 

2.  402 

0.  9772 

2.  518 

306. 00 

0.  4112 

7.  590 

5.  175 

2.  186 

4.  885 

2.  873 

2.  351 

0.  9553 

3.  030 

312. 10 

0.  3589 

7.  647 

5.  885 

1.  937 

5.  040 

2.  850 

2.  311 

0.  9381 

3.  513 

317. 70 

0.  3202 

7.  692 

6.  606 

1.  736 

5.  184 

2.  827 

2.  277 

0.  9233 

4.005 

323. 20 

0.  2888 

7.  739 

7.  443 

1.  550 

5.  333 

2.  804 

2.  258 

0.  9145 

4.  585 

328. 80 

0. 2597 

7.  759 

8.  078 

1.  432 

5.440 

2.  783 

2.  238 

0. 9055 

5.  010 

333.80 

0. 2414 

7.  206 

2.  852 

3.  767 

3.  895 

3.  055 

2.  400 

0. 9829 

1.457 

273.  80 

0.  5873 

7.  264 

3.  080 

3.  517 

4.  025 

3.  029 

2.405 

0. 9842 

1. 607 

278.  80 

0.  5521 

7.  312 

3.  315 

3.  289 

4.  139 

3.  007 

2.  401 

0. 9817 

1.  760 

283.  20 

0.  5195 

7.  395 

3.  714 

2.  968 

4.  308 

2.  981 

2.  397 

0. 9790 

2.  031 

288.  80 

0.  4724 

7.  495 

4.  364 

2.  560 

4.  509 

2.  953 

2.  378 

0. 9703 

2.471 

294. 90 

0.  4107 

7.  590 

5.  196 

2.  178 

4.  693 

2.  931 

2.  343 

0.  9556 

3.  042 

299. 90 

0.  3515 

7.  702 

6.  675 

1.  720 

4.  956 

2.  895 

2.  280 

0.  9283 

4.  060 

308.  25 

0.  2806 

7.  766 

8.  316 

1.  392 

5.  122 

2.  871 

2.  234 

0. 9084 

5.  169 

313.  80 

0.  2287 
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